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ABSTRACT 

The inverse of the fermion matrix squared is used to define a transfer 
matrix for domain-wall fermions. When the domain-wall height M is 
bigger than one, the transfer matrix is complex. Slowly suppressed chiral 
symmetry violations may then arise from all eigenvalues of the transfer 
matrix which are located near the unit circle. Using a variable lattice 
spacing for the fifth coordinate we enforce the strict positivity of the 
transfer matrix for any M. We furthermore propose a modified pseudo- 
fermion action, aimed to decrease the density of close-to-unity eigenvalues 
of the (positive) transfer matrix, at the price of a small renormalization of 
the coupling constant. We explain why these changes may reduce chiral 
symmetry violations in lattice QCD simulations. 



1. Introduction 



Numerical QCD simulations [1-4] using domain-wall fermions [5-8] reveal small, 
but still significant, lattice- artefact violations of cliiral symmetries for A^^ ~ 15, where 
Ng is the extent of the lattice in the fifth direction. It is important to improve our 
understanding of such anomalous effects, and, hopefully, to devise new domain-wall 
actions capable of better suppressing them. Both at tree-level and in perturbation 
theory the anomalous effects decrease exponentially with increasing Ng. This applies 
in particular to the additive radiatively- induced quark mass [P, |], ||. In this paper 
we will thus focus on non-perturbative anomalous contributions to the chiral Ward 
identities introduced in ref. |T^ and first studied numerically in ref. [Q. 

Free domain- wall fermions have a single chiral zero mode on each four- dimensional 
boundary of a five- dimensional lattice when the five-dimensional mass term (or "do- 
main wall height") M is in the range < M < 2. As one moves to the range 
2 < M < 4 this zero mode disappears, and four new zero modes of the opposite 
chirality appear on each boundary. Thus, the range < M < 2 supports a single 
(approximately) massless quark, and the range 2 < M < 4 supports four (approxi- 
mately) massless quarks with "flipped" chiralities. This pattern generalizes to higher 
values of M, until for M > 10 there are no zero modes at all. 

With suitable subtractions to ensure its finiteness (to be discussed in detail in this 
paper) the iV^ — > oo limit of domain-wall fermions can be written down compactly in 
the overlap formalism [|ll|, |l^]. In this formalism the (subtracted) fermion partition 
function is expressed in terms of the overlap of a second-quantized ground state with a 
reference state or, more generally, as the ground-state expectation value of an operator 



representing the boundary conditions in the s-direction [10 



In the "old" overlap formalism [^ , the second-quantized hamiltonian is (minus) 
the logarithm of the transfer matrix that hops domain-wall fermions a single site 
in the s-direction. This formalism has the same massless-quark spectrum as above. 
The "new" overlap formalism [|12] involves a different hamiltonian which arises in a 
continuous s-coordinate limit. The new overlap (and the related four-dimensional, 
non-local, chirally invariant action |jl3|) has the same massless-quark spectrum for 
< M < 2. In the range 2 < M < 4 four massless quarks with flipped chiralities 
appear as before, but the original massless quark remains in the spectrum. As M is 
further increased more massless quarks appear, alongside with all the previous ones 
that keep staying in the spectrum. This different free-field spectrum agrees with the 



pattern of level crossing found in a smooth instanton background [^ 



A more fundamental difference between the two overlap formulae is found in the 
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properties of the transfer matrix. In the new overlap formahsm the transfer matrix 
is strictly positive for all values of M, simply because it is defined from the outset 
as exp{—H), where if is a well-defined local lattice hamiltonian. In the old overlap 
formalism, on the other hand, the transfer matrix is positive only for < M < 1, as 
in the framework of domain-wall fermions from which it is derived. We argue below 
that for M > 1 both positivity and hermiticity of the domain wall's transfer matrix 
are lost, and its eigenvalues are in general complex. 

The aim of this paper is to study what could be the dominant anomalous con- 
tributions to chiral Ward identities, and to suggest methods for suppressing them. 
A key role will be played by a (first quantized) transfer matrix related to the in- 
verse of the fermion matrix squared. For flavor non-singlet chiral symmetries, the 



anomalous term in the lattice Ward identities [|T0| involves the correlations of fermion 
operators on the s = 1 and s = Ng boundaries with fermion operators at s ~ Ns/2. 
Such correlations decrease exponentially if the transfer matrix has no eigenvalue with 
an absolute value close to one. Conversely, if there are many such eigenvalues, one 
expects a much slower (presumably power-law) decrease of anomalous correlations. 

Numerical results ^ ^ suggest an optimal value of M between 1.6 and 1.7 for 
/3 ~ 6. In this range, the transfer matrix of standard do main- wall fermions is complex. 
Slowly decreasing anomalous contributions may then arise from all eigenvalues of the 
transfer matrix which are close to the unit circle. 

The positivity of the transfer matrix can be enforced for any M by restricting 



the range of a^, the lattice spacing for the fifth coordinate ||T5[, a result that could 
have been anticipated in view of the relation between the old and new overlaps. The 
transfer matrix is strictly positive for Mas < 1- (Thus, for M < 2 it is sufficient to 
take e.g. = 0.5.) In this case, only close-to-unity eigenvalues |l^, [16|, |T7 | 

may lead to significant anomalous contributions. We propose a method to reduce 
the density of such eigenvalues in dynamical simulations. The method consists of a 
modification of the pseudo-fermion (Pauli-Villars) part of the action, and it might be 
effective already for modest values of Ng. The main side-effect of the modified action 
is a small renormalization of the coupling constant. 

This paper is organized as follows. Domain-wall fermions with a variable lattice 
spacing for the fifth coordinate are reviewed in Sec. 2. The (first quantized) transfer 
matrix is introduced in Sec. 3. The case of a complex transfer matrix is discussed 
in Sec. 4. The modified pseudo-fermion action is introduced and discussed in Sec. 5. 
Our conclusions are given in Sec. 6, and some technical details are relegated to two 
appendices. 
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2. Domail-Wall Fermions with a variable s-spacing 



Allowing for a variable lattice spacing 05 for the fifth coordinate [T^ , the domain- 
wall fermion matrix is 
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(2.1) 



The four-dimensional lattice spacing a is set to unity. The above matrix structure cor- 
responds to the fifth coordinate s, which we assume to take the values s = 1, 2, . . . , A^^- 
Each entry is a four-dimensional matrix. Pr^l = !(!='= 75) denote chiral projectors, 
and D is the Wilson-Dirac operator 
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M-W C 
-Ct M-W 



(2.2) 



where 



a 



xy 



w. 



xy 



9 5Z \^x+fi,yUx^i Sx-fL,yUy^ 



(2.3) 
(2.4) 



The matrices C and W correspond to the kinetic and Wilson term respectively. We 
also define 

Bxy = (1 - Ma5)6xy + a^Wxy . (2.5) 

[B is proportional to the identity matrix in spinor space. The latter may be either 
two- by-two or four-by-four; the correct meaning can be inferred from the context.) 

Assume momentarily a semi-infinite range for the s-coordinate. For sin(p^) = 0, 
H = 1, . . . , 4, the free-field equation has a right-handed homogeneous solution 



^|JUs■,P,) = = [1 + a,{Wo{p,) - M)Y 



(2.6) 



where the subscript zero denotes free-field quantities, Wo{p^) = I]^i(l ~ cos(p^)), and 
at the corners of the Brillouin zone WqIp^j) takes the values 0,2, ... ,8. The above 
homogeneous solution is a zero mode {i.e. it is normalizable) provided 

-l<l + a,iWoip^)-M)<l. (2.7) 
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We depart from the original domain-wall framework by replacing the constraint on 
the lower bound with the stronger one 

< 1 - Mas • (2.8) 

Since is a positive matrix, this implies 

< 1 + a5(spec(iy) - M) . (2.9) 

The last condition ensures the strict positivity of B. The upper-bound constraint 
in (|2.71) is simply 

Woip^)-M<Q, (2.10) 

which is evidently independent of a^. For < M < 2 there is a single zero mode at 
= 0. If we increase M while decreasing to maintain the constraint ( p.8|) , new 
zero modes will appear at M = 2, 4, 6, 8, while all the zero modes from smaller values 
of M will remain in the spectrum. 

Assuming < M < 2, on a finite lattice the right-handed part of the quark 
field corresponds to the above zero mode, whereas the left-handed part corresponds 
to another zero mode located near the s = Ng boundary. Both at tree level and in 
perturbation theory, for m = there is a small mixing between the two chiral modes 
that vanishes exponentially with A^^ ^, ^, ^]. The parameter m in eq. ( |2.1| ) is 
related to the bare quark mass. An easy way to see this is to invoke the free domain- 
wall hamiltonian. For m = 0, the eigenvalues E{pk) of helicity eigenstates are given 
by = J2k=i sin^(pA:)- For m 7^ we find using first-order perturbation theory that 
E{0) = rrig where 

mg = mMa5(2 - Mas) , (2.11) 
generalizing the result of ref. @|. 



3. Propagators and transfer matrices 

The familiar relation D = 7sZ)^7s, valid for the Wilson-Dirac operator, generalizes 
in the case of domain-wall fermions to 

DF = n^,D^^n^,, (3.1) 

where TZss' = Ss,Ns+i-s' 113 • Eq. (|0| ) implies relations between second-order op- 
erators: DpDp = [TZ'^^DpY = 7lj5{DpDp)7l'y5. We also note that TZ'j^Dp is 
hermitian. For definiteness we will focus on the operator DpDp. One has 

DpDi = ^,B'/^n^,B'/\ (3.2) 
5 



where explicitly 
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(3.3) 



2 + al^,B-'/'DD^^,B-'/' 

Y + PLB-\{ma5)'' - I) 

Y + PnB-\{ma,)^ - 1) 
X ^ = mas • 



(3.4) 
(3.5) 
(3.6) 
(3.7) 



The factors have been introduced for later convenience. Thanks to our insistence 
on the strict positivity of B, these factors are strictly positive too. An interesting 
observation is that 

Y = f + f-\ 

where 

B~' + alB~'/^ CC^ B-'/^ a^B-'/^CB'/^ 
a^B^'^ Ct fi-V2 B 



(3.8) 
(3.9) 



The basic properties of the transfer matrix T follow from 
where 



K 



fi-V2 
asCt 5-1/2 ^1/2 



51/2 



(3.10) 



(3.11) 



which imply that T is strictly positive, bounded, and has detT = 1. 

The first-quantized transfer matrix usually encountered in the context of domain- 
wall fermions is T = KK^ . Evidently, T and T have the same spectrum, and (up to 
normalization) the eigenvectors of T are obtained from those of T by multiplication 
with K. The appearance of T instead of T is due to a technical reason that we explain 
later. 

Eqs. ( p.4|) and ( p.8|) imply that an eigenvalue of T is equal to one iffthe hermitian 
operator 75D has a zero mode. The last condition implies det (D) = 0, an equation 
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which defines a measure- zero subset of the gauge field configuration space. In the rest 
of this section we will assume that det (D) ^ 0, and hence that no eigenvalue of T is 
exactly equal to one. 

We now turn to the construction of the domain- wall propagator Gp = Dp^. Using 
eq. (|3^) one has 

Gp = D^p^,B-'/'G^,B-'/', (3.12) 

where G = Our task is to find an explicit representation for G. We begin by 

writing the spectral decomposition 

f = Y.\vi)\i{vi\ . (3.13) 

i 

For each eigenvalue, we define qi = min(Ai, A~^). We now construct a new matrix 

i 

with the property < spec(Q) < 1. Next we consider the infinite-size matrix 
constructed from the translationally invariant part of Q, extended to the range — oo < 
s,s' < oo. The inverse Goo = is 

Goo{s,s') = ^|t^i)gl'"''/(gi) (wil (3.15) 

i 

= (3.16) 

where f-\q) = q-^ - q > 0. Eqs. (pD and (U) imply Es^oo{s" , s) G^{s, s') = 
for s" 7^ s', while the correct normalization for s" = s' is ensured by the presence 
of f{Q). Since it was constructed using Q (and not T), Goq{s,s') vanishes for large 
\s-s'\. 

Returning to the finite-s case, we now have for 1 < s, s' < A^^ 
G(s, s') = Goo(s, s') + H++{s, s') + H^4s, s') + H+^s, s') + H.+{s, s') , (3.17) 
where 

H__{s,s') = Q^^+^-M__ Q^=+^-^' (3.18) 

The four-dimensional matrices A±± solve a system of linear equations given in Ap- 
pendix A, and have a convergent Ng oo limit. Hermiticity of G implies A ^. = A\__. 
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The physical significance of the above construction is the following. In a fixed 
gauge-field background U = {Ux^} one has Gf{s, s';U) = (ip{s) 'ijj{s')) . Suppose 



'u 

that we project ip{s'), the anti-fermion field on a given s'-slice, onto the (four- 
dimensional) state where Vi is an eigenvector of T. Using eqs. (|3.12|) and 
( p.l7| ), each term in the corresponding projection of <^'?/'(s) -^(s')^^ must have an s'- 
dependence given by either A^** or A^** (Aj is the corresponding eigenvalue). The 
s'-dependence will involve no other eigenvalue of T. A similar statement can be made 
for the ip-Held. To this end, we write Gp = {DpDp)~^Dp and use eq. ( p.l| ) to related 
{DpOp)'^ to {DpDp)^^ . This reproduces a key feature of the second-quantized 
transfer matrix formalism [11, 10 1, but now in the lagrangian formalism, which is 



directly related to the manner one simulates fermionic correlation functions. 

We end this section with two technical comments. In the second-quantized trans- 
fer matrix formalism one can define s-dependent operators via aL,R{s) = T^di jiT'^ , 
where T is the second-quantized version of T (defined by T = exp(a^ logT a), where a 
and are creation and annihilation anti-commuting operators). As noted in ref. ||10[ 
(see eq. (A. 9) there) using the transfer matrix T implies that the operators dji{s) and 
dpls) are identified with the Grassmann variables iPr{s) and iPl{s — 1) respectively. 
In other words, related Grassmann variables and operators do not always have the 
same s-coordinate. On the other hand, if one uses the second- quantized version of 
T, then related Grassmann variables and operators do have the same s-coordinate in 
all cases. This explains the appearance of T instead of T in the expressions for the 
domain-wall propagator. 

Our last comment concerns the new overlap formalism. Dropping all terms of 
order and higher in T one has 

If we now take the limit 05 while keeping the product asiVg = Lg fixed, we have 
T^'ia^ — i> exp(75DLs). This reproduces the result of ref. |jl5|. Recognizing that 



75-D is the new overlap's hamiltonian, we see that both the old and the new overlap 
formulae for QCD can be recovered as suitable limits of domain-wall fermions. 

4. Complex transfer matrix 

In this section we discuss what happens if Ma^ > 1, starting with the free-field 
case. Since [Bq,Cq] = 0, one can lump together B^^^"^ factors, so that only B^^ will 
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occur in the equations that define Qq and Yq {cf. eqs. (|3.2|) and (|3.4|) ). One has 



Yoip,) = B,\p,)[l + + Y: sin^p,)] . (4.1) 

Note that the sign of Yo{pfj) is determined by the sign of Bq{j)^), and that either 
^o(P/i) > 2 or Yq{p^) < —2. The eigenvalues Tq{p^) are all real, and have the same 
sign as Bq{j)^) {cf. eq. (|3.8|)). Therefore, the free-field propagator exhibits sign 
oscillations, but otherwise everything stays pretty much the same as in the Ma^ < 1 
case. 

The situation is different in the interacting theory, where C] ^ 0. When 
Mas > 1, -B has both positive and negative eigenvalues, and B^^"^ has both real and 
imaginary eigenvalues. Consequently, the transfer matrix T is no longer hermitian. 
(Eq. ( |3.1CI| ) still holds if the same definition of 5^^/^ is used in the expressions for K 
and K^^ , cf. eq. ( p.ll| ), even though 5^/^ is no longer hermitian. Of course, K^^ does 



not stand for the hermitian conjugate of K in this case.) 

While det (Dp) is always real due to eq. (p.l|), it is not necessarily positive for 
Mas > 1- In particular det (Dp) may occasionally vanish. (This is not true for 
Mas < 1 and m > 0, see Appendix A.) The propagator exists except on the measure- 
zero subset defined by det (Dp) = 0, and can be constructed using the same technique 
as in Sec. 3. The matrix Q is now defined as follows. Let i/i denote an (in general 
complex) eigenvalue of Y. For yi that does not belong to the closed interval [—2, 2] 
on the real axis, we define qi to be the solution of qi + qj^^ = yi obeying \qi\ < 1. For 
Hi G (—2, 2) the roots obey = 1 while qt ^ q~^ . In this case we arbitrarily pick one 
of the roots. Finally, we disregard gauge-field configurations leading to any yi = ±2, 
hence f{qi) exists. (This amounts to ignoring another measure-zero set; see also the 
last paragraph of Appendix A.) The rest of the construction is the same as before. 

Again, as discussed in Sec. 3, the dependence of the propagator on s and s' is 
governed by the (now complex) eigenvalues and eigenvectors of the transfer matrix 
T. Hence, slowly decreasing anomalous effects may now arise from eigenvalues lying 
anywhere close to the unit circle. 

When Mas > 1, both B~^/^ and T are unbounded. The possibility that T may 
have very large eigenvalues is, however, not worrisome. The latter correspond to 
very small eigenvalues of Q, and so they lead to very short-range correlations in the 
s-direction. 

Of special significance is the chiral Ward identity that governs the pion mass. 
For any as and M, the anomalous term in the pion-mass Ward identity is positive 
when the number of dynamical fiavors is even, as well as in the quenched case (see 



Appendix B of ref. 1|10|)- Thus, the anomalous contributions to this particular Ward 
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identity, coming from all (real or complex) eigenvalues of the transfer matrix, always 
add up. 

Because of eqs. ( |3.4| ) and ( p.8|) , close-to- unity eigenvalues of the non-hermitian 
T should still correspond to approximate zero modes of the hermitian operator 75-D. 
The latter have been extensively studied recently [§, |^. To date, however, no 
information exist on the distribution of eigenvalue in all the rest of the complex 
plane. In particular, it is not known how many eigenvalues are located near the unit 
circle away from the point one on the positive real axis. Therefore it is also not known 
how much of the anomalous effect observed in numerical simulations [1-4] is due to 
(approximate) zero modes of 75-D. The above issues clearly deserve a more detailed 
study. 



5. A modified pseudo-fermion action 

In the rest of this paper we impose the condition Mas < 1, and alongside with 
it the strict positivity of the transfer matrix. For simplicity we also restrict the dis- 
cussion to < M < 2. As discussed in the introduction, the troublesome eigenvalues 
of the transfer matrix are now the ones close to unity. In this section we propose a 
method to reduce the density of close-to-unity eigenvalues in a dynamical simulation 
by modifying the Pseudo-fermion (also known as Pauli-Villars) part of the action. 

Let us first recall why a domain-wall fermion action must be accompanied by a 
pseudo-fermion action. From the point of view of the gauge field, the domain-wall 
action introduce A^^ "flavors" of four-dimensional Dirac fermions. For < M < 2, 
only one Dirac fermion is light (and is identified with a quark field). The other A^s — 1 
Dirac fields have 0(1) masses. If their number was kept fixed, we could simply ignore 
them in the continuum limit. The chiral limit, however, requires A^^ — > 00. If this 
limit is taken at fixed value of the bare coupling (7, the A^s-dependent contribution of 
the heavy "flavors" must be subtracted. 

One can express the domain-wall fermion determinant as 

det {Dp) = fip' X (finite factor) . (5.1) 

The "finite factor", which accounts for the quark field, has a convergent A^, ^ 00 
limit. The bulk term is the (undesirable) contribution of the 0{Ns) massive 
"flavors" . Explicit expressions for both terms were first derived in the transfer matrix 
formulation |Tl], |1^], and more recently by direct manipulations of determinants [p!8| . 
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The latter technique will also be used below. Explicitly, 



= n ' (5-2) 

Ai>l 



where the product is over the greater-than-one eigenvalues of T, cf. eq. (p.9|). As 
expected, the bulk term is independent of the choice of boundary conditions in the 
s-direction [i.e. it is independent of m). One way to cancel the bulk contribution is to 
introduce a five- dimensional boson field (pxs, having the same spin and internal indices 
as the domain-wall fermions, on a lattice whose fifth coordinate ranges only from 1 
to = Ns/2. We will refer to (p^s as the pseudo-fermion field. The pseudo-fermion 
action is 

^unmodified ^ Y: <l>UD,,Dl,Uys4ys' ■ (5.3) 
xs,ys' 

A common choice is Dp{ = Dp{ma^ = 1). Excepting ma^ = —1, in fact, any 0(1) 
value for ma^ will do. When integrating over both fermions and pseudo-fermions 
the bulk terms cancel out, leaving a convergent Ng ^ oo result that has the same 
long-distance behavior in four dimensions as the finite factor in eq. ( |5.1| ). (Another 
option is to use a first-order pseudo-fermion action. When the transfer matrix has 
a gap, the bulk term of the first-order action converges faster to the fermionic bulk 
term This advantage disappears if the gap is small.) 

As mentioned above, for Mas < 1 chiral symmetry violations should arise from 
close-to-unity eigenvalues of the transfer matrix. Actually, in the free-field case the 
spectrum has a gap. This follows from eqs. (|3.4|) and (p.S]), and the fact that 



DoDlip,) = (M - Woip,))' + Esin'(PM) > • (5.4) 

Since the Brillouin zone is compact, the (M-dependent) minimal eigenvalue of DqDq 
is strictly positive. Now, the eigenvalues of DD'^ are gauge-invariant continuous 
functions of the link variables. Since DqDq has a gap in the free-field case, one 
expects that DD"^ too should have a (somewhat smaller) gap if the local plaquette 
action 

V, = Y: Retr (1 - Ux,U,^i,,,Ul^,^^Ul) , (5.5) 

is very small everywhere. In other words, zero modes of DD^ should arise only if 
Vx exceeds an 0(1) constant co(M) > at least for some lattice sites. Hence, if 
one is sufficiently close to the continuum limit and M is not too close to the critical 
points or 2, these zero modes should be suppressed by the plaquette action. In the 
context of domain-wall fermions we can therefore regard these zero modes as lattice 



artefacts. (See refs. [|T^, |T3 for related work. The existence of a co(M) > can be 
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tested numerically. A proof that Co(M) > is desirable, but one should keep in mind 
that rigorous inequalities may under-estimate the value of co(M).) 

For /5 ~ 6, one is probably not close enough to the continuum limit to effec- 
tively suppress the zero modes of DD^^ by the plaquette action alone. The above 
consideration lead us to propose a modified pseudo-fermion action 

-^pf = E <Pls{D^iDlU,ys'<Pys' + E 0L(ciPx)>x. , (5.6) 
xs,ys' XS 

where Ci is a continuous parameter and n is a positive integer. Assuming Cq{M) 
is known, a reasonable choice is Ci ~ Co^(M). The choice of n is discussed later. 
The idea behind the modified pseudo-fermion action is to suppress in a selective way 
those gauge-field configurations supporting close-to-unity eigenvalues. The aim is to 
achieve this suppression in (dynamical-fermion) simulations with not too large A^^ and 
at presently accessible values of (3. As explained below, the main side-effect of the 
modified action is expected to be a small renormalization of the coupling constant. 

As in eq. ( p.2|) , the pseudo-fermion matrix defined by eq. ( |5.6|) can be written as 
75^^/^ fipf 75!?^/^. The explicit expression for Vlpf is the same as the r.h.s. of eq. (|3.3| ) 
except for the replacement 

Y ^ Fpf (ci) = Y + B-"^ icVr B-^/^ , (5.7) 

where V stands for the diagonal matrix d^yVx- (Similar replacements are made in the 

definitions of X++ and X For ma^ = 1 one has Fpf = = XP^_.) 

An expression for det (f^pf) can be written for even N using the general formulae 
of ref. |jl8|, see Appendix B below. The result is 

det (^]pf) = det (R' - PT^^ , (5.8) 



where 



= (5.10) 
Tpr = i " , \ (5.11) 



P= \ I . = (1 - Ql,)-' \ \ T- I . (5.12) 
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As before, the matrix Qpf is defined by Qpf + Qpl = Vpf and the condition < 
spec((5pf) < 1. One has 

PT'P''= I I . (5.13) 

Note that the r.h.s. of the last equation is a function of Ypf only. 

For large N, vanishes whereas Qp/^ grows unboundedly. Like eq. (|5.1| ), one 
can express det (^^pf) as /i^(ci) times a finite factor, where 

/xpf(ci) = det(g;/(ci)) . (5.14) 

(See eq. ( p.7|) for the finite factor.) For A^^ = 2A^, the total bulk factor coming from 
the integration over both fermions and pseudo-fermions is 

(5.15) 

In the unmodified case, Ci = 0, using eq. ( |5.2| ) and det (T) = 1 one has /ipf(O) = 
showing that the bulk factors indeed cancel each other. 

We now discuss how the bulk factor is modified for ci > 0. Since the plaquette 
term is positive, the Ci-derivative of the eigenvalues of Ypf is always positive, and the 
Ci-derivative of the eigenvalues of Qpf is always negative. Therefore /il'Z/ipf (ci) is a 
decreasing function of ci. 

Let us now assume that T (or Q) has an eigenvalue very close to one. For 
M ~ 1.7, the support of the corresponding eigenvector should consist of very few 
lattice sites ITBl. Moreover, as discussed earlier Vx is likely 0(1) on those sites. 




As Ci is varied from zero to its chosen value, we thus expect an 0(1) change in 
the corresponding eigenvalue of Qp{{ci). We will argue below that there should be 
only a small change in the product of all eigenvalues which are not close to unity. 
Consequently, the 0(1) change in what used to be a close-to-unity eigenvalue implies 
an 0(1) reduction in the magnitude of /i|,//Xpf (ci). This, in turn, should suppress the 
Boltzmann weight of the corresponding gauge-field configuration already for modest 
values of Ns {cf. eq. ( |5.15| )). 



Eigenvalues of the transfer matrix not too close to unity are typically not localized, 
and it is plausible that the effect of the modified action on them can be accounted for 
by perturbation theory. It is easy to see that the leading perturbative effect of the 
modified action is to renormalize the coupling constant as 

V<E(^rviAs) . (5.16) 
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Since the free transfer matrix has a gap, the perturbative value of (y'Pxs'Pxsj at any 
finite order, and in any gauge-field background, is regular. (This is true even if the 
associated exact transfer matrix has a unit eigenvalue. In other words, a singularity 
in the propagator cannot develop if we sum the Born series only up to a finite or- 
der.) Now, the expectation value on the r.h.s. of eq. ( ^.16|) involves n or more loops. 
Therefore the resulting change in should be of order NgO^. Since in practice 
a ~ 0.1, if we take for example = 4 this effect may be at the level of 1% or less, 
when Ns is in the range of 10 to 100. Thus, the leading side-effect of the modified 
action seems to be an innocuous, relatively small change in the bare coupling. 

In short, we believe that perturbation theory can be trusted for the collective 
contribution of all eigenvalues, except when there are close-to-unity ones. Since the 
perturbative effect should be small, an 0(1) change in /ipf(ci) should take place only 
when there are close-to-unity eigenvalues, and this change works in the direction of 
suppressing the Boltzmann weight of the corresponding gauge-field configurations. 

The above arguments, while plausible, are heuristic. One question that can 
be settled by an explicit (perturbative) calculation is whether the loop integrals in 
(Vx^'^ 4>is4'xs), while regular, happen to produce large numerical factors. A calcula- 
tion of this expectation value is also necessary in order to be able to compare results 
at different values of Ci, while maintaining a fixed value of the effective bare coupling 
(c/. the r.h.s. of eq. ( |5.1ti| )). 

Last we discuss how the Ng ^ oo limit may be taken with the modified action. 
Evidently, if all other parameter were kept fixed, then in the limit Ng oo the 
perturbative effects induced by the modified pseudo-fermion action would eventually 
run out of control. While the strict A'^, ^ oo limit is not very useful for practical 
purposes, it is legitimate to ask whether, in principle, the modified action has a 
sensible Ng ^ oo limit. 

If we allow both Ng and n to grow, the limit may in fact depend on the ratio of 
these two numbers. As an example, one option is to take the limit n oo before 
the limit Ng ^ oo. Sending n to infinity has the following effect. For ciVx < 1, the 
limit of (ciVx)"' is zero. Hence, the modification vanishes if CiVx < 1 for all x. If, 
on the other hand, CiVx > 1 even for a single site, the norm of (ciVx)" will blow up, 
and along with it some of the eigenvalues of Qp/. In summary, in the limit n oo 
the bulk factor /i|,//ipf (ci) is unchanged if CiVx < 1 for all x, whereas it vanishes 
if ciVx > 1 for any lattice site. Therefore, the limit n oo amounts to imposing 
the constraint ciVx < 1 on the gauge-field configuration space. While this constraint 
should not change the continuum limit, throwing out all configurations with some 
CiVx > 1 in a completely unselective manner could slow down simulations. As argued 
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above, for moderate values of Ng and relatively small values of n, the modified action 
may do a good job in suppressing only those gauge-field configurations supporting 
close-to-unity eigenvalues, leading to a minimal "waste" of configurations. 

6. Conclusions 

In this paper we have studied potential sources for the anomalous term in chiral 
Ward identities. For conventional domain- wall fermions (05 = 1) the transfer matrix 
is complex if M > 1. As discussed in Sec. 4, in the numerical work of refs. [1-4] slowly- 
suppressed chiral symmetry violations could therefore arise not only from eigenvalues 
which are close to one, but in fact from eigenvalues in the vicinity of the entire unit 
circle. 

For 05 = 0.5 and < M < 2 the transfer matrix will be strictly positive. The 
first question that has to be addressed numerically is whether the favorable range for 
a single light quark is still M ~ 1.7, as in the 05 = 1 case 0, §]. The range M ~ 1.7 
seems reasonable also from the point of view of the new overlap formulation [p!6| . 
Since the new overlap involves a continuous s-coordinate limit, it is plausible that the 
best value of M may be rather insensitive to 05. 

For fixed 1 < M < 2, as 05 is decreased from 1 to 0.5 all eigenvalues must flow 
towards the positive real axis. Depending on the flow pattern, the number of close- 
to-one eigenvalues at = 0.5 could be quite different from the original number of 
eigenvalues close to the unit circle. Hence, the transition from 05 = 1 to 05 = 0.5 could 
by itself have a signiflcant effect on the anomalous term. Finally, once the above issues 
are resolved, one can proceed to test whether the modifled pseudo-fermion action is 
useful in further reducing lattice-artefact violations of chiral symmetries at presently 
accessible values of A^^- 

Acknowledgements 

I thank the participants of the RIKEN BNL Research Center workshop on Fermion 
Frontiers in Vector Lattice Gauge Theories, held on May 6-9, 1998, at Brookhaven 
National Laboratory, for extensive discussions that motivated this work. Especially 
I would like to thank the organizers for creating this unique opportunity to exchange 
ideas on the subject. This work is supported in part by the Israel Science Foundation. 



15 



Appendix A. The second-order propagator 



For 2 < s" < A^, - 1, eqs. Q and Q imply 

n{s", s) G^{s, s') - 6s",s' = E ^(^"' ^) H±±{s, s') = 0. (A.l) 

s s 

For s" = 1 and s" = Ng there are boundary effects. The above expressions are not 
zero, and their s'-dependence is given by Q^^' . The linear combination that gives 
the propagator G (eq. (|3.17| )) is determined by requiring the coefficients of Q^^' to 



vanish. By imposing this condition on the Q"^*' part we obtain 

/(Q), (A.2) 
and by imposing this condition on the Q"'^ part we obtain 

/(Q), (A.3) 





- x__g-i 



where 



Note the C has a convergent Ng oo limit. A second-quantized transfer-matrix rep- 
resentation of the finite- fermion determinant is given in ref. [|10| (see in particular 
eqs. (3.1) and (3.10) therein). We decude from it that for Ma^ < 1 and m > 0, 
det (Dp) is (real and) strictly positive for any number of flavors. Therefore Dp has 

an inverse, and the above equations must have a solution which obeys A ^ = A\__. 

We note that Dp^ may exist even if T has an eigenvalue which is exactly equal 
to one. In this possible way to construct Dp^ is to perturb the background 

field so that that eigenvalue will be only approximately equal to one, say to up 0(e). 
The construction of Sec. 3 is now applicable. Dp^ for the initial background field can 
then be found by carefully removing the perturbation, keeping track of the leading 
terms in e. 



Appendix B. The second-order determinant 



For even A^, we can use the formulae for the determinant of a general tridiag- 



onal matrix, derived in the appendix of ref. [18|, to write down an expression for 



det (fipf(ci)), cf. eq. (|5.8|) . (Here the block entries of the tridiagonal matrix have spin 
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indices ranging from one to four. With the replacement Ng, a similar formula 

holds for det (fl), cf. eq. (|3.3|).) In the notation of ref. [0, the present application is 
defined by 



a 

aN/2 




j = l,...,iV/2-l 



N/2 - 1 



(B.l) 
(B.2) 
(B.3) 
(B.4) 



With these definitions, eq. (A. 10) of ref. [T^ reads 



det {a 



det 



-1 / -1 -1 



N/2 



(B.5) 



Substituting the above explicit expressions we arrive at eq. ( |5.8|) . Note that {—a~^l3 
is equal to the r.h.s. of eq. ( p.l3| ). The bulk factor of det (^^pf) can be separated as 
follows. We first rewrite eq. (|5.8|) as 



(B.6) 



det = det (l - Q%) det (S' - T^, 5+) , 
where = P~^R^. Using eq. ( ^.11| ) it is now easy to check that 



det {VL 



pfJ 



det (q/) det (l - Q% 



X det 



Q^f 



s- 



(B.7) 



The bulk factor det (Qp/^) appears explicitly in the above equation. The other terms 
are by definition the finite factor. For Ma^ < 1 the entries of f2pf are bounded, and 
det (f^pf) is finite for any finite A^. This is more easily seen using eqs. ( |5.8|) and (|5.13| ). 
It can be checked that if an eigenvalue of Q approaches one, then the restriction of 



the r.h.s. of eq. ( 5.13 ) to this eigenvalue has a finite limit. Moreover, the restriction 
of PT^ P^^ to this eigenvalue grows only linearly (not exponentially) with A^. 
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